Effect of strain on stripe phases in the Quantum Hall regime 
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u ■ Abstract 

Spontaneous breaking of rotational symmetry and preferential orientation of stripe phases in 

a • 

H ! the quantum Hall regime has attracted considerable experimental and theoretical effort over the 

'^ ' last decade. We demonstrate experimentally and theoretically that the direction of high and low 

^ ■ resistance of the two-dimensional (2D) hole gas in the quantum Hall regime can be controlled by 

an external strain. Depending on the sign of the in-plane shear strain, the Hartree-Fock energy 

K^ ' of holes or electrons is minimized when the charge density wave (CDW) is oriented along [110] or 

^^ ' [110] directions. We suggest that shear strains due to internal electric fields in the growth direction 

cn 

in 

o 

o 
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are responsible for the observed orientation of CDW in pristine electron and hole samples. 



Interplay between kinetic energy and electron-electron interactions in two-dimensional 
electron gases in magnetic fields leads to a rich variety of possible ground states, ranging 
from incompressible Lauglilin liquids, the Wigner crystal, charge density waves (CDW) to 
exotic non-Abelian anyonic states. The possibility of the formation of a CDW state had been 
suggested l| even before the discovery of the quantum Hall effect, and later it was predicted 
that a CDW should be the ground state for partially occupied high Landau levels[2|, yj. 
Experimentally, anisotropic magnetoresistance (AMR) and re-entrant QH phases have been 
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observed in 2D electron[^, |5| and hole[6|, |7| gases. The majority of experiments have been 
conducted on samples grown on (001) GaAs. Unexpectedly, the CDW was found to be 
consistently oriented along [110] crystallographic direction in these samples, a surprising 
fact considering isotropic nature of the wave functions on the high symmetry (001) surface. 
Search for the physical origin of the broken symmetry and the observed preferential ori- 
entation of stripe phases has been actively pursued experimentally and theoretically over 
the past decade. Reduced symmetry of the interface was suggested [8|] as a factor which 
introduces an anisotropy of the effective massjst] or of the cyclotron motion |lO|. However, 
single-particle effects associated with these anisotropics seem unlikely to be responsible for 
the lar ge m agnitude and strong temperature dependence of the resistance IfJ . Later work 



showed 



121] that the precise symmetry of the 2D gas confining potential is also unimportant. 



and micron-scale surface roughness does not correlate with the stripe orientation. There 
have been theoretical suggestions 13l-ll5| that anisotropic correction to electron-electron in- 
teractions arising from elastic and piezoelectric effects can be responsible for the resistance 
anisotropy. While the free energy of the CDW is minimized in the vicinity of [110] and [110] 
directions, those theories cannot explain why these directions are inequivalent. Progress has 
been made in understanding the effect of the in-plane magnetic field, which has been shown 



to influence orientation of stripes 



16 



17] . the effect being explained by the field- induced 
anisotropy of the exchange potential 18|, [l9|. However, naturally existing preference for 
[110] orientation of the CDW in purely perpendicular field, the same for electron and hole 
samples, remained unresolved. 

In this work we show both experimentally and theoretically that strain breaks the cylin- 
drical symmetry of electron-electron interactions in magnetic field and results in a preferred 
orientation of the CDW. We show experimentally that externally applied shear strain can 
enhance or reduce anisotropy of the resistance and switch low and high resistance axes. Our 
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FIG. 1: Magnetoresistance is plotted as a function of B for the current aligned with [110] (red) 
and [110] (blue) directions in strained and unstrained samples. In a) no external strain is applied, 
(b) thermally-induced tensile strain is along [110], and along [110], £p = 0. Inset shows sample 
schematic, red and blue arrows show current, and green arrows show strain. 

theory shows that depending on the sign of the in-plane shear strain, the Hartree-Fock en- 
ergy is minimized when the CDW is oriented along [110] or [110] directions. We suggest that 
shear strains due to internal electric fields in the growth direction caused by mid gap Fermi 
level pinning at the sample surface are responsible for the observed preferred orientation of 
the CDW in pristine electron and hole samples. Finally, by applying uniaxial strain we are 
able to induces a stripe phase even at a filling factor u = 5/2, with a CDW winning over 
other QH states. 

Samples were fabricated in the van der Pauw geometry from carbon doped GaAs quantum 
well heterostructure grown on (001) GaAs [7|, |20||. From the low field Shubnikov de-Haas 
oscillations, the hole density is 2.25 x 10^^ cm~^, and the mobility 0.8 x 10^ cm^/Vs is 
determined at the base temperature lOmK. Some samples were thinned to 150 /im and 
glued on a multilayer PZT (lead zirconate titanate) ceramic actuator with samples' [110] 
or [110] crystallographic axis aligned with the polarization axis of the PZT. Application of 
voltage Vp to the actuator induces in-plane shear strain in the sample Sp/Vp = 2.8 x 10~^ 
V~^ and small uniform bi-axial strain. The total shear strain e = eth + £p also includes 
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FIG. 2: Strain dependence of the anisotropic magnetoresistance. (a-d)Magnetoresistance in the 
vicinity of u = 7/2 and 5/2 for /||[110] and /||[110] as a function of voltage on the piezoelectric 
actuator Vp. In (e-f) magnetoresistance at i^ = 7/2 and 5/2 is extracted. On the top axis, Vp is 
converted to the piezo-induced shear strain in the heterostructure; actual strain includes thermally- 
induced offset. 

a residual strain eth due to anisotropic thermal coefficient of the actuator, which depends 
on the Vp during cooldown. To insure that voltage on the actuator does not induce charge 
modulation in the attached sample we insert a thin metal foil between the sample and the 
PZT. The foil was also used as a back gate to adjust 2D gas density which has a weak 
dependence on strain due to difference in piezoelectric coefficients of GaAs and AlGaAs 21 1 
(density changes by 3% for the maximum voltage span on the PZT). 

Magnetoresistance for pristine (i.e., not attached to an actuator) sample is shown in 
Fig. [1^. States at z/ = 7/2 and 11/2 are highly anisotropic with low resistance direction 
along [110], while states at z/ = 5/2, 9/2 and 13/2 are almost isotropic, consistent with 
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the previous study j7|. In Fig. [T]d,c similar traces are shown for large e = e^/j ^ +10 ^ 

{sp = 0). For Eth > the anisotropy is enhanced compared to 



and e 



^th 



< 
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the unstrained sample, with states at z/ = 5/2, 9/2 and 13/2 becoming anisotropic and 
resistance for /||[110] approaching zero for half-filled Landau levels. For Sth < the low and 
high resistance axes are switched. Here strain also leads to strong anisotropy at z/ = 5/2 



with high resistance axis along [110] direction. 

Residual strains eth in Fig. [T|d,c are larger than the in situ adjustable strain Sp. In Fig. [2]we 
analyze AMR for a sample cooled with Vp = —150 V, aiming for eth ~ 0. Magnetoresistance 
as a function of Sp is plotted near u = 5/2 and 7/2. At V^ < magnetoresistance is highly 
anisotropic. For the high resistance direction, the resistance -R7||[iio] strongly depends on Sp 
and decreases by a factor of 50 (4.4) at u = 5/2 (7/2) as Vp is varied from -300 V to 300 V. 
i?/||[iio] increases only 1.7 (1.3) times. At Vp > 100 V, the resistance at z/ = 5/2 is isotropic, 
with no maxima for either current direction, as in unstrained sample. Thus, in the range 
of small strains, the CDW is not a ground state at z/ = 5/2, consistent with observations 
in unstrained samples. From the data we conclude that e < within the adjustable range 



of Vp, because i?/||[iio] > -R/||[iio] at z/ = 7/2. Continuous evolution of Ri 



consistent with continuous change reported in in-plane magnetic field 17] 



110] and -R/||[iio] is 



Having presented experimental results of strain on the resistance in the QH regime, 



we now develop the Hartree-Fock theory of CDW by extending theoretical work 22| to 
anisotropic 2D systems. Strain results in sizable modification of the 2D single-particle spec- 
tra, and of the many-body wavefunctions. For either electrons and holes, the 2D Hamiltonian 
in magnetic field H = V x A is 

^ ^ {Px - f ^x) _^ {Py - -c^x) ,^s 

2mx 2my 

where a:||[110], ?/||[110] are the principal axes of the reciprocal mass tensor, with m^^ = 
{m^^ + my^)/2 and yU^^ = {m'^^ — m^ ^)/2 being isotropic and anisotropic parts. For holes 
in a III-V material quantized along the (001) direction irT^ = —Oji + 72 + a73)/?7io, where 



7i;72;73 are negative constants defining the bulk hole spectra 23|, |2J| and the numerical 



coefficient a is defined by these constants 25|] . For the anisotropic part induced by the shear 



strain e = exx = —Syy our result is fi^^ = 'y3\/3de/[y2{Tch/a)'^], where d is the deformation 
potential 23|, and a is the quantum well width (we assume an infinite rectangular QW). 
For electrons, m is the 3D effective electron mass rric, while for the anisotropic part our 
result, obtained in the third-order perturbation theory, is ^^/2/i = —P'^de/\^E'^, where P 



is the Kane band coupling parameter 26| and Eg is the band gap. Note that the sign of the 



strain-induced term for electrons is opposite compared to that for holes. 



To find the single particle wavefunctions we define new coordinates, x' = xJuix/m' , 



y' = yJ^/'^' 1 where m' = .JrfT^m^. In these coordinates, 1-120 = 2^(p'~f ^0^ is isotropic, 



and the wavefunctions are the usual LLs wavefunctions with degeneracy in the guiding 
center coordinate X' = k' J"^, where / = {hc/eHY^'^ is the magnetic length. According to the 



physical picture of Aleiner and Glazman [22| , the electron-electron interactions in a system 
with partially filled high topmost LL can be considered as interactions via the Coulomb 
potential with an effective dielectric constant defined by the electrons of the fully filled LLs. 
In the deformed coordinates, for A^ ^ rj^ ^ 1, where A^ is the LL index, r^ = (vrna^)"^, 
n is the 2D carrier density, and as = h /me^ is the Bohr radius, the low energy physics of 
the 2D electron liquid in weak magnetic field is thus described by the effective Hamiltonian 
^e// ~ Tir ^q' P(^')'^(^')p(~^')' where L^, Ly give the size of the sample, 

^ q = — 2 

-oK(q')v/(a^ + K^)^ 
is the Fourier component of the renormalized electron-electron interaction potential, kq is the 
background dielectric constant, and ^(q') is the effective dielectric constant. An expression 



for K valid at l/(v2iV+T/) < Q' < kp, where kp is the Fermi wavevector, is /€(q') = 1 + 
2/gaB[27|. Note that this ^(q') is isotropic in the original coordinates but anisotropic in the 
deformed ones. Finally, p(q') is the Fourier-component of charge density of the partially filled 
LL N, p(q') = Ex'aN{q')e-^''^^'''"^y''/'^alax'-,'^i^, where a^(q') = L%{^) expi-^), 
L^(x) is the Laguerre polynomial, and 0^1 (o-x') is the creation (annihilation) operator for 
a hole with guiding center at X' in the topmost LL. Defining the CDW order parameter 
^(q') = TIT ^x' <^x'+q' i'^ /2'^x'-q' P /2i wc obtaiu the Hartree-Fock energy as 

Ehf = ,r^ E [«^(q') - «e.(q')]A(q')A(-qO, (3) 

27r/ UN q,_,o 

where < z/jv < 1 is the filling of the topmost LL. The Fourier transforms of the Hartree 
and exchange potentials are, respectively, uh{(() = v{c()[aN{c()Y and 

UeM) = ^j:MQ')^^p[^i'i-QWy + QW^)]- (4) 

J^x-I^y Q/ 

Note that the two potentials are related by a Fourier transform but the arguments in the 
Fourier transform of the Hartree potential are transposed relative to the arguments of Uex- 







While this transposition is not important in the isotropic case [2|, taking it into account 
here is crucial for finding the preferred orientation of the CDW modulation. 



a 



Applying an analysis similar to that in [2|, we see that when the L%{x) in aAr(q') is zero, 
Mj:/(q') is also zero, so that Ehf < 0, as given by the exchange contribution. The system 
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FIG. 3: a) The ratio of isotropic and anisotropic parts of the hole exchange potential (Eq. (H) for 
three values of shear strain e. b) Self-consistent calculations of the band profile and internal electric 
fields in the studied wafer, c) Modeling of a HIGFET structure from f2S ]. 



becomes unstable with respect to the formation of the CDW with a wavevector close to q'. 
However, in contrast to the isotropic system, in which the direction of the CDW is chosen 
spontaneously, in the presence of strain the smallest q', which corresponds to the direction 
with the largest mass, gives the largest value of the exchange and the lowest HF energy. 
In Fig. [3^, we plot the ratio of the anisotropic and isotropic parts of the hole exchange 
potential at d = —bAeV, 71 = —6.8, 72 = —2.1, 73 = —2.9, and a = 0.4. Note that here the 
anisotropy reaches 3% for strains of 10~^. The CDW near the half filling of the A^-th LL 
naturally results in stripes with alternating z/^r = and u^ = 1, which, in turn, translates 
into low resistance direction along the stripes and high resistance direction perpendicular to 
the stripes. Thus, in the presence oi e ^ 0, the CDW has a preferential direction defined by 
the sign of the strain, consistent with the experimental results. While the theory is valid, 
strictly speaking, in high Landau levels it appears to describe experiments even at A^ = 1, 
the lowest LL for which L%{x) has a zero. 

The preceding analysis suggests that internal strain may be responsible for the observed 
orientation of stripes in pristine samples, where no external strain is applied. GaAs is a 
piezoelectric material and any electric field in the z direction results in an in-plane shear 
strain e = di^^Ez-, where c/14 = -2.7- 10~^° cm/V. A calculated band diagram for our samples 
is shown in Fig. [3b. Inside the QW Ez < 2-10^ V/cm and results in strain too small to orient 
the stripes. E^ on both sides of the QW, caused by doping, is also small and, in our samples, 
odd in z. However, in all GaAs samples there is a large field near the surface of the wafer 



due to the pinning of the Fermi energy near mid- gap. This surface charge-induced field is 
~ —10^ V/cm and the corresponding strain is e ~ 3 ■ 10^^. If transmitted to the QW region, 
this strain has the correct sign and magnitude to explain the observed orientation of stripes 
in pristine samples. To show that transmission of strain does indeed occur, we consider the 
minimal model in which this effect is present. The free energy of the model is given by J-" = 
J-'ei+2/3i5j;e, where the elastic free energy is J-'ei = |A [{dxUy + dyUxY + {dzUxY + [dzUyY], Ux 
and Uy denote displacements, e = ^{dxUy + dyUx) is the shear strain, a;,|/, z now correspond 
to the [100], [010], [001] directions, (3 is the piezoelectric constant, and A is the elastic 
one. In realistic devices, Ez{x,y,z) is nonuniform in {x,y) plane due to charge fluctuations 
on the surface and in the doping layer. For illustration, we consider a cylindrical sample 
of radius L and height d with cylindrically symmetric Ez{r,z) = J2nEn{z)Jo{(lnf), where 
r = {x'^ + y'^y^'^, Jo is the Bessel function, and Qn are quantized by the condition jQ^qnL) = 0. 
The solution for the strain is 

^(^' ^) = I: E ^^E^i^)MQnr) + £5(r, z), (5) 

ZA „ o^ q^ 

where e^ is localized near r = L. If Ez{r,z) varies with r at some characteristic scale R, a 
typical Qn is of order 1/R. Then, according to Eq. |5l the corresponding component of the 
strain propagates largely undiminished over distances z ~ i? from the region where En{z) is 
large (sample surface). As a result, macroscopic regions with sizable strain exist throughout 
the QW region. 

The presence of internal strain also explains orientation of stripes found in 2D electron 
gases. A typical band diagram of an electron sample is similar to that shown in Fig. [Sjo, but 
is inverted relative to the Fermi level with both the surface electric field and the shear strain 
changing sign. However, the anisotropic term in Eq. [1] also has opposite sign for electrons 
and holes, so that the sign of the anisotropic term in electron and hole exchange is the same. 
Thus, for both holes and electrons, a surface field will orient the CDW along [110], as seen in 
experiments. With our model, we can also explain the reorientation of stripes as a function 
of density observed in a HIGFET (Heterojunction-Insulated Gate Field Effect Transistor) 

. At low gate voltages (low densities), shear strain will be dominated by the surface 



field, see Fig. [3b. At large gate voltages (2 V corresponds to 3 ■ 10^^ cm~^), the electric field 
across the AlGaAs barrier becomes large enough to change the sign of the strain in the 2D 
gas region, thus reorienting the stripes. 
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In summary, we have shown experimentally that resistance anisotropy in the QH regime 
can be controlled by external strain. Theoretically, we have traced this effect to a strain- 
induced anisotropy of the exchange interaction and a competition between the internal and 
external strain, which defines the preferred direction of CDW. In general, any factor that 
brings in a crystallographic anisotropy of the effective mass gives rise to a crystallographic 
anisotropy of the Hartree-Fock energy of the CDW state (this is the case, for example, 
in hole gases grown on low-symmetry (311) GaAs [6]). Our analysis suggests that, for 
heterostructures grown in the high-symmetry (001) surfaces, piezoelectricity due to surface 
electric fields is the largest source of such anisotropics. We underscore that, although the 
anisotropy of electron Hartree-Fock energy is two orders of magnitude smaller than that for 
holes, it must still choose a preferential direction for the CDW of guiding centers. Therefore, 
the preferential direction of the resistance anisotropy in pristine (001) samples appears to 
be universally dictated by internal strain. 
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